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Abstract 



A formulation for transport in an inhomogeneous, interacting electron gas is described. 
Electronic current is induced by a constraint condition imposed as a vector Lagrange 
multiplier. Constrained minimization of the total energy functional on the manifold of an 
arbitrary constant current leads to a many-electron Schrodinger equation with a complex, 
momentum-dependent potential. Constant current Hartree-Fock and Kohn-Sham approx- 
imations are formulated within the method and application to transport for quantum wires 
is developed. No appeal is made to near equilibrium conditions or other approximations 
allowing development of a general ab initio electronic transport formulation. 



In recent years, there has been considerable advances towards the development of accurate 
theoretical models to deal with electron transport through atomic and molecular wires P[. 
i ' This activity has been largely spurred on by development of experimental techniques to form 
■ atomic scale electrical contacts such as scanning tunneling microscopy [0] and mechanically 
O ■ controllable break junctions ||. Measurement of current- voltage characteristics have been 
H \ recently performed on a single benzene- 1,4-dithiolate molecule [f|], buckminsterfullerene ||, 
.£h I individual atoms ||, and DNA strands |7]]. The mechanisms for charge transport through these 
^ ' low- dimensional structures remain largely unexplored and unexplained. Several approaches 
o5 ' have been developed on the basis of the extended Hiickel || and density functional theory 
within a Lippmann-Schwinger formalism || or nonequilibrium Green's functions |fL0[j . In these 
approaches, the system is partitioned into three parts: molecular device (or wire) and two 
electron reservoirs. 

The interaction between contacts and wire represents the kernel for the Lippmann-Schwinger 
equation, and serves generally to define the non-interacting system for Green's function ap- 
proaches. Recent studies indicate the exact nature of the molecule-contact bonding are critical 
in predicting the correct order of magnitude for the currents with applied bias revealing that 
accurate electronic structure methods must be incorporated into quantum transport methods 
for accurate predictions and analysis of current- voltage relationships. 

In this letter, we present a variational approach to charge transport in correlated electron 
systems. We begin by introducing a many-electron Schrodinger equation which provides exact 
wavefunctions on the manifold of a given current distribution. The many-electron Schrodinger 
equation is developed without any restrictions to the treatment of electron-electron correla- 
tions, whether as injected current or part of the molecular device. We next present specific 
approximations corresponding to the Hartree-Fock and configuration interaction treatments of 
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the many-electron problem. Next, the Kohn-Sham approximation to density functional theory 
is extended to the inhomogeneous electron gas with a fixed current at zero external magnetic 
field. 

We begin by developing a general many-body theory for an inhomogeneous electron gas 
with arbitrary and static current density. For an inhomogeneous, interacting iV-electron gas in 
an external potential v(r), the ground state energy is given by the following energy functional 
(atomic units are employed throughout the paper unless otherwise mentioned): 

E[p,T] = = / h (r)p(r,r')\ r=r ,dr + l [ f r(n, r 2 ;n,r 2 )dndr 2 (1) 

J Z J \T\ — r 2 | 

where ^o( r i) is the one-body electron operator which includes the electronic kinetic energy and 
the one-body electron interactions with external potentials. We have introduced the one- and 
two-electron reduced density matrices: 

N 



p(r;r') = N /*(r, r 2 , . . . , r N )**(r', r 2 , . . . , r N ) f[ dr, (2) 

J i=2 

N 

r(r, r'; r, r') = N(N - 1) j tt(r, r', . . . , t n )V*(t, r', . . . , r N ) J] dr u (3) 



where \1/ is the many-electron wave-function 

The fundamental quantity in our approach is the electron current density 

j(r) = ^[V r -V r ,]p(r,r')| r=r ,. (4) 

The current density is constrained to have specified values at prescribed points in position space, 
or more generally a functional G[j](r) is required to be zero at defined positions r. The task is 
to minimize the energy functional E[p, F] subject to the constraint C[j](r) = 0. Altyhough the 
constraint formulated in this manner is nonholonomic, it may be included into a variational 
functional via pointwise vector Lagrange multipliers a(r). The constraint is explicitly achieved 
by introduction of an auxiliary functional 

Q[p, T, j, a] = E[p, T] - /i (/ dvp(v) -Nj+ A[j, a] (5) 

The first two terms are standard in the variational formulation of quantum many-body theory 
with the first giving the total energy and the second is introduced to constrain particle number 
(or likewise, to introduce the orthonormality constraint for the wavefunction). The third term 
has been introduced to impose the current constraint. The functional A[j,a] depends linearly 
on a(r) and is designed to insure the following condition is maintained: 

Variation of Q[p, T,j,a] with respect to the Lagrange multiplier a provides an additional con- 
straint equation. One can view the a(r) as an infinite set of Lagrange multipliers, each associ- 
ated with C[j](r) = at a point r. The variation of A[j,a] over the current density j(r) yields 
a vector field: 
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Extrema of the auxiliary functional Q[p, T,j,a] corresponds to the following many-electron 
Schrodinger equation 



E [hofo) + j. [V„ A(r,)] + } + g j-l-j j $ = M (8) 

where we have defined the energy E = N/j,. The anti-commutator term, 

[V,A(r)] + = VA(r) + 2A(r)V, (9) 

gives rise to an additional single-electron imaginary potential arising directly from the constraint 
on the current density This additional constraint potential forces the many-electron wave 
function \I/ to be irremovably complex and enforces the required current density distribution at 
extrema. 

A physical interpretation of the vector field A(r) follows if we re-write the eq.(|8]) in the 
following form 

A(r,)) 2 - iA(r,)A(r,) + vfo)} + £ — L-f ] * = (10) 

A } i>j \ r i T j\ ) 

The current constraint has introduced terms equivalent to a vector potential of an external 
magnetic field |p| , however, the vector field is a functional of the \I/ and has to be dtermined 




self-consistently from the Eq.[TD] in our case. 

The constrained many-electron Schrodinger equation is not yet in a form allowing for a 
solution to be found as the the Lagrange multiplier a(r) and the vector field A(r), which is 
a functional of a(r), for a specified current are not given. An expression for a(r) can not be 
found until an explicit form of constraint functional A[j,a] is first introduced. The motivation 
for the present study is to desire to develop the basic many-body formalism for the calculation 
of a direct current I through a quantum wire. To achieve this aim, we specify the constraint 
on the current density distribution in the following form: 

dydzj x (v) = I x (11) 

Within this description, net current flow is aligned along the x— axis. This is a simple geometric 
arrangement to specify current flow for quantum wire and can be easily extended to more 
complex topologies. For the case specified, the functional A[j,a] takes the following form: 

A[j,a] = J dxa x (x)( J dydzj x {v) - I x ) (12) 

The vector field induced by the presence of the constraint on the current density has spatial 
dependence and nonzero projection only along the x-axis: 

A (r) = ^M = M*),0,0) (13) 

Therefore, for a quantum wire with a direct current the vector field A(r) coincides with the vec- 
tor Lagrange multiplier a(r). Generally, if the constraint functional A[j, a] is a linear functional 
of j(r), then the vector field A(r) and the Lagrange multiplier at each point r are equivalent. 
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We next solve for the Lagrange multiplier a x (x) and obtain 



(*) = aT {/ dydz I n + - 2 ^(^)} 



(14) 



where we have introduced the quantity 

p x (x) = J dydzp(r). (15) 

Suppose now that the many-electron wave-function \l/ is approximated as Slater determinant 
of N orthonormal single-electron orbitals ^(r). The use of the Slater determinant enables us 
to write the electron current density (|4]) as the sum of electron current through all occupied 
orbitals 

1 occ 

j( r ) = ^E(^( r )V^(r) -^(r)V^(r)) (16) 

Likewise, as for the derivation of the exact many-electron Schrodinger equation, the current 
distribution is fixed via G[j](r) = 0. The constraint is included into Hartree-Fock minimization 
via an additional Lagrange multiplier A[j, a] = A[ip, if)*, a]: 

Q[M*,a] = E[^*]+A[^rM (17) 

Minimization of eq.flTTD subject to the orthogonalization condition J dri/;*(r)i/;j(r) = 5ij results 
in a Hartree-Fock equations with fixed current density distribution 

Hhf + \ % [V, A(r)] + ) Vi(r) = E^r) , (18) 

where Hhf is the Hartree-Fock operator which contains the electron kinetic energy, external 
potential, Hartree potential and Fock exchange operator, A is again an additional vector field 
induced by the current constraint. The configuration interaction (CI) representation follows 
directly if we substitute for the many electron wave function \1/ a linear superposition of Slater 
determinants or spin coupled determinants. The CI expansion is substituted into eq. || and 
variation with respect to the expansion coefficients yields a modified matrix eigenvalue problem, 
whereby the interaction matrix elements are supplanted with the additional one-body potential 
terms arising from the constraint potential ^ [V, A(r)] + . 

We next turn to the task of incorporating the current constraint into the density functional 
theory (DFT) ||12|| . We describe an extension of the Kohn-Sham formulation of DFT [13 



to systems with a non-zero current density distribution. For an inhomogeneous, interacting 
electron gas in external potential v(r), the ground state energy is given by the Hohenberg- 
Kohn energy functional 



E HK [p\ = T [p\ + J dvv{x)p{v) + P -^l + E xc [p] (19) 

The first term, To, is the kinetic energy functional of iV noninteracting electrons with given 
density and current density distribution. The second and third terms, are the interaction energy 
with external potential and the electrostatic interactions. The last term, E xc [p], is the exchange 
and correlation energy functional. 
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We follow closely to the derivation of the many-electron Scrodinger equation described in 
the beginning of this letter. The distribution for the current density vector is fixed through 
a functional G[j](r) = at certain points r. The constraint is included into DFT variational 
determination of the charge density via a Lagrange multiplier a(r): 

Q[p, j, a] = E HK [p] - p (J dvp{v) -N) + A[j, a] (20) 

Following the Kohn and Sham approach fl3| , we introduce a reference fermion system with 



orthonormal single-particle orbitals ipi(r) and occupition numbers rii to reproduce the charge 
and current densities: 

p(r)=$>^*(r>Mr) (21) 

i 

j(r) = ^^(^WV«r) - V*(*)V#(r)) (22) 
In analogy to the Kohn-Sham approach, we consider variation of the auxiliary functional 



Q[p, j, a] eq. 20 with respect to the single-electron orbitals ipi(r) to obtain the following single- 



electron self-consistent equations: 



~ (-zV + A(r)) 2 - J A(r)A(r) + v(r) + J drj^r + p xc [p](r) 



^i(r) = E^r), (23) 



the desired result. 

In order to get an indication of how our approach works, we solve the Kohn-Sham equa- 
tion with fixed current ( p3| ) for the case of a uniform electron gas. Consider an electron gas 
in a cubical box of volume V = L 3 , throughout which a background positive charge is uni- 
formly spread out rendering the system neutral (jellium model). A direct current is applied 
along the x-axis I = (I x , 0, 0). The classical Coulomb electron-electron interaction, the Hartree 
term, is completely compensated by the electron interaction with positive charge density and 
by the electrostatic energy of the positive background. The vector field is constant for the 
translationaly invariant uniform electron gas yielding the following Kohn-Sham equation for 
the model: 

(--A + fe + -Av)^(r)=E^ t (r) (24) 

with specification of the constraint taken in the form given by eq. (|Tl"D . A dispersion relation 
which is coupled to the constraint equation for the current density may be obtained by inserting 
a plane-wave solution ^(r) = -^exp(ikr) into eq2. (p4] , |TTD , with the result 



k 2 

E k = — + p xc + A x k x (25) 

2 occ 

T E^ = / - ( 26 ) 
L k x 



From the system of equations eqs.( [25| - |26| ) the exact expression for the vector field A can be 
deduced: 

A = -J> ( 27 ) 
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with electron density p = N/L. Given the vector field, we can write the current density 
dependent dispersion relation of the uniform electron gas with applied direct current I x : 

k 2 I 

E ]l = — + p xc - —k x (28) 
z p 

The presence of the constrained current lifts the degeneracy of the k and —k states and results 
in the energy gap between the electron moving in the direction of net current flow, i.e. with 
positive k x , and electrons moving in the opposite direction. 

We finally demonstrate application of the method with a simple numerical example. We 
consider a one-dimensional system with current injection into a fixed, external potential. Trial 
states are taken to be pure plane wave states and these states are then allowed to relax to 
account for the presence of an external potential, taken to be a simple Gaussian form. In fig. 
1, the solutions to the Schrodinger equation are shown with the current fixed by the constraint 
potential, i.e. no boundary conditions are imposed other than through the constraint condition. 
Within the figure, the changes in real and imaginary components of the wave function relative to 
a pure plane wave state are plotted. The constraint re-arranges the solution to the Schrodinger 
equation allowing for the changes due to the current encountering the external potential, while 
maintaining a constant current fixed at the value specified by the constraint condition. 

In this letter, we have given a variational formulation of quantum electronic transport. 
There is no resort to any imposed conditions other than to constrain the current distribution. 
Thus the formulation is equally valid for highly correlated systems and for all nonequilibrium 
current regimes. The formulation has been specified for general many-body theory and cast into 
a form suitable for the Hartree-Fock and configuration interaction methods. We next showed 
how to introduce the constraints into density functional theory. The method as such is quite 
general and applicable to all common approaches to electronic structure theory. 

For the case of one-dimensional transport in a two-terminal system, we have introduced 
a physical constraint which allows for an explicit determination of the Lagrange multipliers 
required to fix the current flow. With the Lagrange multipliers in hand, it is possible to deter- 
mine the complex potential needed for solution of the many-body wavefunction, or electronic 
density, satisfying the physical boundary conditions introduced by a constant current flow in 
and out of a region. The formulation of the problem has been chosen for greatest compatibility 
with the methods in common use for electronic structure theory determination, and results in 
a completely variational formulation of the quantum electronic transport problem. 

Acknowledgment This work has been funded by the European Union through the Infor- 
mation Society's Technology (1ST) programme, within the Future and Emerging Technologies 
Advance Research Initiative's NANOTCAD project (IST-1999-10828). 
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Figure caption 

figure 1- Wave function shifts from an attractive Gaussian potential V ext = V exp(— (x — fi) 2 /a 2 ) 
are shown as the difference in the wavefunction \l> = + it; relative to a plane wave state \1> in 
box normalized units. All other quantities expressed in atomic units. Incident energy is E = 
2.0 a.u., well depth Vo=-0.1 a.u., well breadth a = 0.047T. Shift in the real component — 0o 
are displayed with the heavy lines, shift in the imaginary component £ — £o are displayed with 
the lighter lines. 
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